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1.0 INTRODUCTION 


The REEDS 'owputer code Is a tool useful for the estimation of certain 
rocket exhaust efflu.nt concentrations and dosages and their distributions near 
the earth's surface following a rocket launch event. Output from REEDS has been 
useful In producing near 'real-time air quality and environmental assessments of 
the effects of certain potentially harmful effluents, namely HCL, AI 2 O 3 , CO, and 
NO. REEDS Is the last In an evolutionary series of codes which treats the 
atmosphere within the first few kilometers of the earth's surface as a single 
layer; thus the name REEDS - Rocket Exhaust Effluent Diffusion Single layer 
model . 

REEDS was developed for the NASA tropospheric environmental program by 
the NASA/Marshall Space Flight Center (MSFC) Environmental Effects Task Team 
(EETT) member. Science Applications, Inc. NASA leadership was provided by the 
EETT leader, Dr. J. Briscoe Stephens. As a result of continuing research, the 
function of REEDS Is being largely replaced by the more comprehensive multi-layer 
diffusion program, REEDM. Therefore, this document should valuably assist the 
process of evaluating the progress made by the application of new techniques to 
the assessment problem. 

REEDS Is designed to operate on a Hewlett Packard 1000 series computer 
system. Research, development and testing took place on the MSFC system at the 
Space Sciences Laboratory (SSL), followed by Implementation and subsequent 
operational monitoring support and bio-medical studies via the system at the 
Kennedy Space Center (KSC). A modest effort was made to keep these systems 
compatible for modeling software operations purposes, even though their diverse 
purposes obviously gave the systems different capabilities, peripheral hardware, 
and software. 

An ^siiportant feature of the HPIOOO series computers used In the REEDS 
code Is the Program Scheduling feature. The feature permits one self-contained 
program to Initiate another completely self-contained program, and then return 
operation to the first. This Is somewhat analogous to the FORTRAN feature of 
calling SUBROUTINES. Thus, a program which completely fills the available 
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program memory may execute another, which may execute another, etc., greatly 
expanding the capabilities of the computer system. Thus, the main executive 
module schedules major modeling modules, which. In schedule other modules, 
facilitating a logical division of research and mode .4 effort. The resulting 
overall structure In the REEOS application Is Illustrated In Figure I-l. 

REEDS was originally developed for interactive graphical operation 
through an SAI Plasmascope^*. The requirement for compatibility with the KSC 
installation necessitated the development of a graphics handler permitting 
operation with or without the Plasmascope^. This handler uses several assembly- 
level progranned subroutines, which constitute the only non-FORTRAN software in 
the package, and which permit complete operation via the HP Cathode Ray Tube 
(CRT) terminals as well as the Plasmascope^. 

While the J1EEDS code Is set up primarily for operation on Titan, Delta, 
and Space Shuttle Vehicles at their respective KSC/Cape Canaveral Air Force 
Station (CGAFS) sites. It has a capability for far nw-e general applications. 
The vehicle data, primarily source rate and altitude histories, are Input in 
parametric form. Simple replacement of these and the site coordinates allows the 
user to apply the code to a far more general class of problems. This Is 
especially true since the code can handle both the continuous and instantaneous 
sources. 

Further documentation regarding structure and modeling will be found 
in Section 2 of this document. Operational guidance Is contained In Section 3, 
and program and input data file listings will be found In Section 4. 


*Tra demark of an SAI Technology Division Product. 
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SECTION 2.0 

REEDS MODEL JESCRIPTION AND BASIC STRUCTURE 


The overall structure of the REEDS code is broken down into blocks 
which largely mirror the phenomenological categorization of the effluent 
diffusion processes. There is first a "Thermodynamic Phase" in which a rocket 
exhaust cloud is generated, transported and expanded until approximate 
thermodynamic equilibrium with its environment is reached (Figure 2-1). The 
RCLDS module (cloud rise module) contains most of the modeling for this phase. 
RCLDS calls RMTPS, which illustrates the local meteorology on a plotter for the 
user, showing approximate geometric relationships with respect to the effective 
atmospheric boundary layer interface. Then, thermal equilibrium having been 
reached, the turbulent diffusion process is allowed to proceed (Figure 2-2) in 
the "Kinetic Phase". The RCONS module contains most of the modeling for this 
process, and plots the concentration and dosage as a function of range along the 
cloud centerline. RCONS also calls RISOS, which plots concentration and dosage 
contours (called "ISOPLETHS") on a map of the launch vicinity. Of course, 
control, meteorological, and configuration input is provided by the main REEDS 
executive .Module, as is sequence control. 

To make input more flexible, use is made of "question" input files. 
For each module for which it is necessary, user input prompts are stored in an 
indexed file, from which a module may draw in order to prompt the user for a 
particular input. Changes and/or additions may then be made as necessary to 
these files. Listings of these files are shown in Section 4. 

The structured diagram of Figure 2-3 illustrates the gross operational 
sequence of the entire code. This will be elaborated upon in the following 
subsections. 

2.1 MODULE DESCRIPTION 

Each module will now be roughly described in narrative fashion, 
seasoned with the appropriate mathematical expressions employed in the modules. 
As much as possible, the narrative "walk" through each module will follow the 
commenting in the code. Some commenting not directly related to calculational 
procedures may be excluded from the description. In addition to making use of 
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FIGUnF 2 2. ‘LLUSTRATION OFTHE *'KIMF1 1C” PHARE MODELED IN REEDS 



Figure 2-x structured diagram op the reeds code showing Important and commented 
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FIGURE 2-3. (Continutd) STRUCTURED DIAGRAM OF REEDS CODE 
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I the coimienting in the code as "landmarks", use will be made of the structured 

diagram of Figure 2>3 for guidance in module description. Each major 
calculatiunal and structural event or user input option will be noted as a 
paragraph heading throughout these descriptions. 

Common Data . The REEDS code makes use of bulk or "blank" common to 
simplify the use of subroutines. A penalty was paid for this in the form of a 
slight loss of modularity and some pain and agony due to vulnerability to "side 
effects" (indirectly and inadvertently changing a variable during a 
modification /upgrade). The penalty was lessened by structuring the conroon into 
testable "pseudo" blocks (not to be confused with block or labeled common), and 
by the use of descriptive comments for each significant "pseudo" block. This 
data is shared by the module and some of its subroutines, and is written to a 
"cownon" disk file. The common disk file is created by the program with the name 
"■REEDS" when the first attempt is made to write the common data. It is finally 
destroyed (purged) by the program upon successful termination of the program. If 
an abnormal termination leaves the file unpurged, it is, in principle, written 
over when it is detected to pre-exist; however, some strange results have 
occasionally been, produced when the file pre-existed, prompting prudent users to 
be sure that it was purged before execution. Thus, following a brief statement 
of purpose, a date for last modification and any necessary loading instructions, 
a large block of coimon data (316 words) appears in each module. Just before 
another module is called this data is written on the disk file. Then, as soon as 
the new module is in core, the disk file is read, reloading the common. When the 
new module concludes, or prepares to call still another module, the common data 
is written again. A brief description of the common variables (Figure 2*4) 
appears in Table 2-1. 

The REEDS "Question" Files . Several REEDS modules have their own 
question files, within which are stored informative messages and/or user 
prompts. These are indexed sequential coded (and therefore listable) files, 
program access to which is accomplished by calling subroutine "DREAD" (for Disk 
READ . The arguments to DREAD are: 

1) A six (6) character, three-word file name 
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TAttt 2 * 1 common data ottcnimof^ 


*muoo- 

COMMON 

•LOCK 

AHOQHAM LASIU 
eaUIVALf NCI NAMI 
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DfSCRIPTION 
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KAOOIt 
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VAAA<1)/0C1 
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TOTAL SOURCE OUTPUT RATE (t/Mc). NORMAL 
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■UMN 
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TOTAL SOURCE OUTFUT RATE {%/mol ASfiORMAL* 




SINGLE ENGINE SURN 


VMK(3>/QC3 


TOTAL SOURCE OUTFUT RATE ASNORMAL* 




EXALOtlON/tURN 


VrAR<4>/QT1 

°C 1 

TOTAL SOURCE STRENGTH (p), NORMAL 


VPAA<6)/QT2 


1 
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^2 

SINGLE ENGINE 


VI*AR<f)/QT3 


TOTAL SOURCE STRENGTH (t), ABNORMAL* 
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•lit ROCKET RISE POWER LAW COEPPICIENT 
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b 

• 2 nd ROCKET RISE POWER LAW COERPICIENT 


VRAR(9>/CC 

c 

•3rd ROCKET RISE POWER LAW COEFFICIENT 


VPARdOWHEATN 


EFFECT HEAT RELEASE (cal/p), NORMAL BURN 

j 
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EFFECTIVE HEAT RELEASE (cM/«L ABNORMAL* 




SINGLE ENGINE BURN 


VPAR(12)/HCATA 
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TAM.1 2-1. ComWNMd 


muoo 

PnOQRAM LIVtL 



COMMON 

CQUIVACKNCi 



BLOCK 

NAMf 

SYMBOL 

oiscRirriON 

omoN 

PAKAMBTERS 

YP08 

— 

CURRENT TERMINAL DISPLAY CURSOR ORDINATE 


YES 

— 

LOGICAL ‘YET* CODE (• 1 MY) 


ISNOFO 

— 

MfTSOROLOQtCAC DATA TYPE INDICATOR 




0 - SOUNDING 

1 - FORECAST 


IVHICL 

— 

ROCKET VEHICLE INDEX 




(1, 2, 3, 4, 8 POR SPACE SHUTTUE. TITAN, DELTA, 
DELTA 3914, AND MINUTEMAN, RESPECTIVELY) 


(RUN 
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PROORAM EXECUTION MODE (OPERATIONAL PRODUCTION) 


(CALC 
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CALCULATION FLAG 


(PLACE 

— 

METEOROLOGICAL DATA LOCATION CODE 


* (TIME2 

— 

METEOROLOGICAL DATA TIME ZONE CODE 


(GOTMP 

— 

FLAG TO INDICATE WHETHER OR NOT A MAP ORiGIN HAS BEEN 
SET 


NUMRUN 

— 

TOTAL NUMBER OF RUNS (MET DATA FILE TO PROCESS) 


LU 

— 

INTERACTIVE TERMINAL DEVICE CODE 


NUM 

— 

RUN SEQUENCE NUMBER 


MONLY 

— 

METEOROLOGICAL PLOT SWITCH ^ 
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1 



1 - YES 

0 - NO 




USED IF ONLY PURPOSE OF RUN IS WET PLOTS 


MPORM 

— 

METEOROLOGICAL PLOT FORM SWITCH 



i 

1 - YES 

0 • NO 


LUPRNT 

— 

PRINTING LOGICAL DEVICE KEY 




6 - PRINTER OUTPUT, ELSE NO PRINTER OUTPUT 


(TOU 
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TOWER DATA SWITCH 
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1 



0 • NO TOWER DATA 

1 - TOWER DATA 
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IPU 
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METEOROLOGICAL PLOTTING DEVICE CODE 
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IPL2 
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CENTERLINE PLOTTING DEVICE CODE 


IPL3 
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tSOPLETH PLOTTING DEVICE CODE 
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TABLE 2-1. ContiniMd 


PSEUDO 

COMMON 

BLOCK 

PROORAM LEVEL 
EQUIVALENCE 
NAME 

SYMBOL 

DESCRimON 

LAUNCH 

TIME 

LTIME 

— 

LAUNCH TIME PLAQ (TRUE IF LAUNCH TIME READ; OTHERWISE. 
FALSE) 


LTIM 

— 

LAUNCH TIME OF DAY (0001; 24001 


LDAY 

— 

LAUNCH DAY OF MONTH (01; 31 1 


LMONTH 

— 

2 WORD ARRAY FOR LAUNCH MONTH OF YEAR (t. 9 ., 2HJA, 1 HNI 


LYEAR 

— 

LAUNCH YEAR (t. 9 ., 1981) 


LAUNTD 

— 

ARRAY OF CHARACTERS COMPRISINQ THE LAUNCH DATE AS 
READ FROM THE TERMINAL 

CLOUD 

STBALT 

H 

STABILIZATION HEIGHT (METERS) 

STABILIZATION 



AND LAYER 

STBHT 

X 

C/I 

m 

EFFECTIVE STABILIZATION HEIGHT (METERS) 

PARAMETERS 



STBAZD 


AZIMUTH FROM PAO OF STABILIZATION (DEGREES) | 


STBRNG 

Rt 

RANGE FROM PAO OF STAPILIZATION (METERS) 


STBTIM 


TIME FROM LAUNCH OF STABILIZATION (SECONDS) 


CLORAO 

Rc 

RADIUS OF CLOUD (METERS) 


RISTIM' 


30-WORO ARRAY OP TIMES FOR PASSAGE OF CLOUDS THROUGH 
EACH MET. LAYER 


BOTLAY 


HEIGHT OF BOTTOM OF LAYER CONTAINING STABILIZATION 
ALTITUDE 


TPLAY 

1 

Zt 

TOP OF LAYER CONTAINING STABILIZATION 


CALHT 

z 

CALCULATION HEIGHT (Zg X Z < Zj) j 


LAYTOP 

— 

INDEX OF LAYER TOP ARRAY ELEMENT CONTAINING I 

STABLIZATION | 


LAYBOT 

— 

INDEX OF LAYER BOTTOM ARRAY ELEMENT CONTAINING I 

STABILIZATION ! 

j 


REFLEC 

T 

SURFACE REFLECTION COEFFICIENT i 

I 

0 • TOTAL REFLECTION i 

1 - TOTAL ABSORPTION ^ 


(PT2 


LEAST SQUARE FIT FLAG FOR POTENTIAL TEIviPERATURE 1 

GRADIENT CALCULATION | 

1 - NO FIT 

0 - USE THE FIT 


SAI.127M 
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TAtUM. ComtfHKd 


^fUOO 

PROQAAM CIVIC 



COMMON 

lOUIVACINCI 



BLOCK 

NAMI 

SYMtOL 

OESCRfPTION 

MCT DATA 

ACT 

Z 

ANMAY OP CAYin-TOP ACTITUOIS IN THI ItilT PROPICI 


lOIK 

9 

AANAY OP WIND OIRICTtONS OP THI CAYIRS OP THI m2T 

PROPtCI 


S^CEO 

u 

ARRAY OP WIND SPEEDS OP THE LAYERS OP THE MET PROFILE 


TEM^ 

T 

ARRAY OP TIMPIRATURIS OP THI CAYIRS OP THI lv.IT 

PROPICI 


PRESS 

P 

ARRAY OP PRESSURES OP THE LAYERS OP THE MET PROFILE 


PTEMP 

PT 

ARRAY OP POTINTIAC TIMPIRATURIS OP THI CAYIRS OP 
THf MCr PROPICI 


SUROEN 


SURPACI DENSITY OP THI MIT PROPICI 

1 

PILE 

— 

PtCI NAME ARRAY (6 CHARACTIRS. 3 WORDS) OP THE MET! 
PICE 

CALCULATION 

ITIME ' 

...M 

CACCUCATION TIM(( OP DAY (0001; 2400) 

TIME 

lOAY 


CACCUCATIQN DAY Qt'' MONTH (01:31) 


1 MONTH 

— 

2 WORO ARRAY FOR CACCUCATION MONTH OP YEAR 
(•.«., 1981) 

! SOUNCiNO 

ISTIME 

— 

SOUNDING TIME OF OAY 

TIME 

♦ 

ISOAY 


SOUNDING DAY OF MONTH (01:31) 

1 

i 

1 

i 

ISMONTH 

— 

2-WORO ARRAY FOR SOUNDING MONTH OF YEAR (•.«., 
2HFE, 1 HB) 

i 

ISYEAR 

— 

SOUNDING YEAR (•.g., 1981) 
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i PSCUOO 

PROGRAM LEVEL 



, COMMON 

EQUIVALENCE 



1 tlOCK 

NAME 

SYMtOC 

DESCRIPTION 

! OtPFUSION 1 

SIGXO 


STAIi UZIO CLOUO ALONQ WIND STANDARD OtVIATION | 

; COiPFtCIINTS 



1 

1 

[ 

Stax 

oX 

ALONO-AXIS t CLOUO EXTENT | 

i- 

SIGAP 


WIND AZIMUTH STANDARD OSVIATtOK AT STAtiUZATION 


XOIST 

X 

ALONG AXIS DISTANCE 

i 

YOlST 

Y 

ORR'AXIS DISTANCE 


EXP2 


VERTICAL DOSAGE TERM 

i 

SQStGZ 

ZoZ* 

TWICE THE SQUARE Of THE 1- VERTICAL CLOUD EXTENT 


SI6Y 

OY 

Of P AXIS l-CLOUO EXTENT 

' 

SIGAL 

^AL 

LAVSRcD-AVERAQfiO WIND AZIMUTH STANDARD DEVIATION 

1 

CLOLNG 

1 

c» 

k. 

CLOUD LENGTH IM)RTERS) 

1 

AVSPO 

\ 

u 

SPEED OF CLOUD AVERAGED OVER MET LAYERS IMETERS/ | 

SEC) 


1 AVOIR 

A 

DIRECTION OF CLOUD AVERAGED OVER *iiET LAYERS ! 

(DEGREES) ! 

1 


SI6AZ 

>^A2 

WIND AZIMUTH UNCERTAINTY (DEGREES) \ 

i CONCENTRATION 
AND DOSAGE 

j CONMAX 

1 


MAXIMUM PEAK CONCENTRATION (ffMI 

! VALUES AND 

CONCPK 

\p 

PEAK CONCENTRATION (ffMI 

: RANGES 


r 



DOSMAX 

°M 

MAXIMUM PEAK DOSAGE (fPM-SEC) j 

1 

1 

DUSPK 

°P 

PEAK DOSAGE iPPM-SEC) | 

i 

! 

RNGSMC 

%c 

RANGE OF MAXIMUM CONCENTRATION (FROM PAD) 1 

1 

i 

1 

RNGSMD 

”mo 

RANGE OF MAXIMUM DOSAGE (FROM PAD) 1 

. - - - - . — 1 
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2) The desired line nuii^er (Index + 1) 

3) A sixty-four (64) character (32 word) data buffer to 

receive the ASCII data. 

The first line of each question file Is an eighty (80) character sequence of the 
numbers 

1234567890123456789012. . . 

useful for preparation/editing of the question file Itself (I.e., aids In 
columnlatlon of chf'racters), so the actual line number Is always one (1) greater 
than the Index of useful data to be retrieved from the question file. Section 
4.3 shows the REEDS question file listing. Special CRT display enhancement 
characters are often embedded as part of the character strings. They are 
unprintable with the selected print option used In generating the listings of 
Section 4, and do not show up even as blank spaces. This gives an erroneous 
column Indication for characters which do show up. However, the special 
Character Display feature of all HP CRT's permits display of the enhancement 
character symbols, facilitating editing of the lines of these files. Once read 
by OREAD, these lines are printed on the desired device by subroutine "CHAR", 
which performs the necessary terminal -dependent oerations required to display 
the line read. 

Cotmiunication of Data To (Or From) Output Files . When it is 

desirable to communicate program numeric data to (encoding) an output file, or 
from (decoding) that file, it must be first internally converted to (or from) an 
ASCII character string. Encoding is accomplished with a call to the system 
subroutine "CODE" and a subsequent specialized "READ" operation. Decoding 
requires the use of "WRITE", instead of "READ". The optional argument to CODE is 
the number of characters desired to be encoded (or decoded). Many painful 

agonies were avoided by explicitly defining and using this argument. Also great 
debugging horrors were generated by inadvertently putting some other statement 
between the call to "CODE" and its subsequent "READ" (or "WRITE"), or otherwise 
failing to resolve a call to "CODE" with the necessary "READ" (or "WRITE") call. 
(These are mentioned here in the event users use this document as a reference in 
making any changes to REEDS.) The arguments to "READ" (or "WRITE") are normally 
the standard output device variable and the FORMAT code. However, when "READ" 
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(or "WRITE") 1s used to resolve a call to "CODE" to encode (non-ASCII) or decode 
(ASCII) data, the first argument Is the array variable name used to receive the 
encoded (or send the decoded) data. The size of this array must be con^iatlble 
with the FORMAT statement (the 2nd argument) used to accomplish the "READ" (or 
"WRITE"). Once encoded, the ASCII string Is written to an output, whenever it is 
desired, in the normal fashion, i.e., with an "A" FORMAT formatted "WRITE" call. 
If the data were decoded (via the "WRITE" call) instead, the resulting data 
(FORMAT permitting) could be arithmetically cr otherwise processed, and must be 
re-encoded if the ASCII string is needed. 

2.1.1 The REEDS Executive Module 

The REEDS executive module 

Initializes some global operational, terminal 
control, and calculational data 

Reads additional control and parameteric data 

- . Reads and prepares the meteorological data 

Calculates the potential temperature profile 

Controls subsequent scheduling of modules. 

The bulk of the module consists of interactive option selection and related 
parameter initialization. The main program (file "REEDS") is 1157 source lines 
long at this printing, and the subroutine packages required (files "REEDT", 
"LINE", "DUNE", "CALIX" and "CTLB2") constitute another 1555 source code lin«s. 
Refer to Figure 2-3 as necessary to follow the program structure. 

User Terminal Identification Setup . The type of terminal (CRT, 
Plasmascope) is first discerned by the program from system tables internal to the 
HPIOOO operating system. Some terminal parameters are subsequently initialized 
~ automatically and transparent to the user. The header is transferred from the 
REEDS executive module question file (called "»REEDS") to the terminal display. 

Cloud Rise Fit or Virtual Potential Temperature Lapse Rate Fit 
Option . The desired fit to the Lapse Rate of Virtual Potential Temperature is 
then produced for the cloud rise model as a result of a user option. The choice 
is between a deterministic method or a least squares method suggested by 
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EETT leader Or. J. B. Stephens. The latter inethod produces consistently lower 
altitudes than the former for cloud rise stabilization and has thus been termed 
the **low altitude fit". This choice remains as a user option pending documented 
validation of one method or the other. This option was installed after 
Plasmascope operation was largely phased out. Thus, the normal terminal 
independent input sequence is not used, rather, a direct CRT terminal input Is 
used. 

Number of Runs and Input Meteorological Data File Names . The REEDS 
program was designed to process multiple meteorological data files if that is 
desired. The procedure is to set up a series of data files which have in conmon 
the first four (4) characters of their names, and a relative sequence number, 
i.e., 01, 02, . . . 99, for their last 2 characters. When prompted to do so, the 
user enters: 


XX, FFFFYY 

for the number of runs, data file name root, and code YY, respectively, where 

♦ 

XX » Number of runs desired 1 01 ; 99 1 

FFFF » The first four (4) characters of the series 

of meteorological data file names 

YY * The file of the series it is desired to start 
with, c=-g., OATAOl 

The program has a default data file name root of "DATA" which may be invoked by 
either assenting to the use of the default (which also assumes there will be only 
1 run of DATAOl), or by entering blanks for FFFF ( which causes any entry for YY 
to be ignored, presetting its value to 1, thus causing the series OATAOl 
OATAXX to be executed). One run of an arbitrary file name is accomplished by 
declining the default, and entering 

XX, FFFFYY 
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1 




I 


where XX ■ 01, and FFFFYY Is the six (6) character file nante desired, e.g., 

01, TESTFZ 

Geographical Location and Time Zone Codes . As a result of a standard 
NASA naming convention the geographical location, and therefore the time zone 
code, can be embedded within the meteorological data file name. These are the 
common variable ITIfiEZ ind IPLACE, respectively. The following illustrates the 
implications of certain non-standard values for FFFF: 

FFFF ■ VANO (Vandenberg AFB), IPLACE » 1, ITIMEZ » 

IHP 

FFFF • KSC (1965 data tape), IPLACE ■ 2, ITIMEZ * 

IHE 

FFFF » KSC, IPLACE - 3, ITIMEZ » IHE 

FFFF « SPEC, Special location, IPLACE and ITIMEZ « 

. 0, and actual subsequently input according to 

the format: values 

:XX,YYZ,AAAA 

where 

XX » Day of month 101 ; 311 

YYY » Month IJan ; Decl, and 
AAAA = Time 1 0001 ; 24001 . 

Any other input is assumed to be the standard ICSC meteorological data file with 
ITIMEZ * IHE, and IPLACE ■ 2. The meteorological data file can only be an actual 
sounding for FFFF * SPEC, whereas for ether inputs sounding and/or TOWER data may 
be used. 

0 

Research. Operational, or Production Modes . As of the time this 
document was being generated the only validated mode of these alternatives was 
the operational mode — which is also the default and most commonly used mode. 
While the additional user inputs inherent to the research mode will be mentioned, 
their demonstration in Section 3 will have to await a subsequent revision of this 
document. Nevertheless, the appropriate non-default selection is made by 
entering a single character R, 0, or P, respectively. This results in the 
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setting of cornnon variable IRUN to I, 2, or 3, respectively. Default selection 
Is, of course IRUN ■ 2. 

Atmospheric Mixing Layer Top and Wind Azimuth Uncertainty - For 
Production Runs . For production runs, the Input of the top of the atmospheric 
mixing layer, TOPLAY, (In meters) and the standard deviation of the wind azimuth 
uncertainty, SIGAZ, (In degrees) must be Input from the keyboard. The format, 
upon prompt Is 

XXXX.XX 

for TOPLAY, with possible range I0;999999.0l meters, followed by 

XX 

for SIGAZ, with possible range (0;991 degrees. For the other operational modes 
(Research, Operational) the values of TOPLAY and SIGAZ are estimated and then 
optionally changed by the user. This will be elaborated upon when it Is 
encountered In the RCLOS cloud rise module (Section 2.1.2) 

Meteorological Forms Switch . The REEDs program plots some of the 
meteorological data (called a "Met Plot") and a schematic of the stabilized 
exhaust cloud. The form for this plot can also be prepared by the program on the 
plotter, but the process is very time consuming. Therefore, it Is better done as 
infrequently as is necessary and high-speed reproduction devices used to 
generate a supply of the forms. For the purpose of generating this form, the 
REEDS program may be run to the point of this input for the sole purpose of 
generating the form. The run may then be terminated and subsequently restarted 
once a supply of forms is available. Thus, when prompted for "Met Plot Forms?", 
the user may respond "Y" for 'Yes' or "N" for 'No', the latter of which is the 
default response. This sets the comnon (OPTION pseudo block) parameter MFORM to 
1 or 0, respectively. 

Meteorological Plots . The meteorological plot may also be skipped if 
it is desirable. When prompted, a "Y" or "N" indicates 'Yes' or 'No', 
respectively. This sets the common OPTION pseudo-block) parameter MONLY to 1 or 
0, respectively. The first setting is the default setting. 
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Line Printer Output . The line printer output may be omitted If 
desired. When prompted a “Y" or "N“ Indicates 'Yes' or 'No', respectively. This 
sets the common (OPTION pseudo-block) parameter to 6 or 64 ^q respectively. The 
first setting Is the default setting, and Its value of 6 Indicates the standard 
printing device code. The value of 64 ^q for the alternate setting formerly 
cleared the parameter (largest value was 77g • 63 ^q). Now, however, the device 
code range has been extended above lOOg, making 64^0 a legal device code. 
However, no such code Is active on the MSFC or KSC machines, so no change (e.g., 
64 jQ — 128 j^q) has yet been made. 

Launch Time and Date . The current computer system time and date are 
loaded Into the launch time and date variables, LTIM, SOT, LDAY, LMONTH, LYEAR 
via the call to the subroutine GETDT, which reads the system clock. The user is 
prompted for this information, and if the response to the prompt is "N" for 'No', 
the system tim^i/date is rejected as the launch time, and opportunity 1s given for 
the above variables to be read accordingly to 

TTTbbbbbODbMMMbYYYY 

where TTTT » the time of day (0001; 24001 , 

b or bbbbb » blank space place holders Indicating the number of user- 

input blank spaces 

DD * the day of the month 101;311 

MMM « the month IJAN;DEC1 

YYYY - the year (e.g., 1980) 

If the response to the original prompt is "Y" for Yes or blank, the system time is 
accepted for the default launch time/date, and the time/date variables are left 
as initialized. 

Launch Vehicle . The user Is then prompted for the launch vehicle 
desired. The options are "D" for Delta, "T" for Titan, and "S" for Shuttle. This 
sets the common (OPTION pseudo-block) variable IVHl'CL to 0, 1, or 2, 

respectively. 


Meteorological Data Input . If the IPLACE parame' .r, set by the 
meteorological data file name read above, is 2, then the forecast/sounding data 
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Is assumed to reside on tape and is read from device code 8 via subroutine KSC65 
(explained below). Otherwise the appropriate disk file name Is opened. Whether 
the data represents an actual sounding or a forecast Is Indicated by a code read 
from the file header. If a sounding Is Indicated, the data type parameter, 
ISNOFO, (an OPTION pseudo«b1ock common variable) Is set to 0; If It Is a 
forecast. It Is set to 1. As the header and/or data lines are read, they are 
printed on the printer device. The following Indicates the sequence or data 
lines read. 

(1) test code, data type and location header Information 

(2) time of data, ascent and data header Information 

(3) actual data for each altitude: 
altitude (feet) 

wind direction (degrees from north) 
wind speed (knots), 
temperature (°C), 

Dew Point (°C), 
pressure (millibars) 
relative humidity (percent) 

Absolute humidity (gms/M^) (if present on file) 

Air density (gms/M ) (if present on file) 
other data (if present) (Ignored by program) 

(4) Terminator (NNNN) indicating end of data. (See Section 4.4) 

The program is designed to ignore all non-essential data words and/or data lines. 
The altitude and wind velocity are then converted to other units; they are 
converted as follows: 

altitude — m 


wind speed — m/sec 


All data lines with altitudes less than 10000 feet are then sorted by ascending 
order of altitude. The potential temperature is then calculated according to 



1 


PTEMP^ - (TEMP^ + 273.6)/lOOO 

y J 

where PTEMP, TEMP, and PRESS are the COMMON variables (MET data pseudo-block) 
potential temperature, temperature, and pressure, respectively. The sorted data 
1s then written on the output device and the meteorological data file closed. 

Met File Reads Via Subroutine KSC65 . The subroutine KSC65 handles 
MET file reads from tape (device code 8). The sequence of data Is Identical to 
that of the disk file read. The subroutine was designed to read large data 
bases, especially the 1965 KSC meteorolocical data tapes (thus the subroutine 
name). The arguments to KSC65 facilitate selection of specific data sets from 
the tape. They are as follows: 

IBUF ■ Eighty (80) character (40 word) integer buffer array to hold 
ASCII data lines read from device 8 pending decoding. 

lALT » Altitude data array to be passed back to the main program for 
' eventual conversion to the analogous common array. 

DPTEMP ■ dew point temperature array to be passed back to the main 
program 

IFOFFl ■ (KSC65 variable IWANT) ■ ordinal data set indicator for 
selection among data sets of the same date/time. 

lEOF » End of File (EOF) flag set if EOF encountered on device 8 
(set by ( ulpment Table Check Three system EXEC call) 

IT IMG ■ Sounding time desired 

IDAYG « Sounding day desired 

IMONG * Sounding month desired 

ISNDG ■ data type desired (sounding/forecast) 

RUNNUM * run numbers desired 

Creation of the Common Data Disk File . The common data, though not 
completely defined yet, is written to the disk file ''»REEDS" via a call to the 
subroutine "RWDIS." The two arguments to RWDIS are the six (6) character common 
data disk file name (3 integer words, 2 characters each), and the operation 
integer code (1 * write). 
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RWOIS contains the 316 word coimon block as It exists In the main 
program. In addition, space Is provided for the data control block array 
(locally called ODCB) necessary to control the disk read (transparent to user - 
only provision of the space necessary, I.e., 144 Integer words or more for the 
control array). The data Is transferred to/from disk via a buffer array. While 
the common data transferred Is 316 words, disk storage Is only by Integer word 
equivalents. Thus, because of the real data varlables/arrays, the Integer word 
equivalent of the 316 words Is 557 words. The buffer array Is thus 557 words. 

The ISIZE array, which Is the parameter controlling the file creation operation, 
must have Its two Integer words set to 30 and 557, respectively. The first word 
Is the number of blocks the 557 words are written In. 

The file Is first opened. If It Is not detected to pre-exist, the file 
Is created with the GREAT system call. For a RWOIS operation code of 1, the file 
Is written; if 0, the file is read. The file is then closed. 

Transfer to the Cloud Rise Module . As can be seen In Figure 2-3, 
transfer to the Cloud Rise Module RCLDS then occurs. The transfer is 
accomplished via the system EXEC call, which schedules the program whose name Is 
contained In a six (6) character (3 word) argument array (RCLOR). The program 
must have an ID segment in CORE, an operation accomplished via the File Manager 
RP command which must be executed before the running sequence begins (see File 
Manager or Terminal Users manual). 

Program Termination . At the conclusion of RCLDS, transfer is 
returned to REEDS and the user is prompted to decide for program termination ("Y" 
or default response) or for the processing of more data C'N'* for No termination). 
If the latter response is chosen, REEDS begins again. Otherwise, the CRT screen 
1s cleared, and the conwon disk file is purged. Finally, all program ID segments 
in the REEDS series are removed from core ("OFF'ed" -- the reverse of the "RP" 
operation — see File Manager or Terminal User's manual). The program STOP is 
executed and REEDS Is terminated. 

2.1.2 The RCLDS Cloud Rise Module 

The RCLDS Cloud Rise Module models the thermodynamic phase of the REEDS 
code. In so doing, the meteorology, energy, and potential temperature factors 
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dre calculated. Then the range from the pad of the cloud center, the cloud 
transport direction, stabilization height, and surface mixing layer depth are 
calculated. These are necessary for the formation of one of the more significant 
calculations -- the wind azimuth uncertainty. This parameter Is an Important 
driver of the turbulent diffusion process. Once the surface mixing layer height 
Is estimated and optimally refined, the RMTPS module Is called to plot the 
pertinent geometric relationships and display the input meteorology and 
calculated potential temperature. The program is then ready for the kinetic 
phase of the calculation. 

Data File Reads . The prompt for the surface layer height Is 
transferred from the RCLDS question file, ''?RCLD$'* (Figure 4), to a holding array 
via the DREAD subroutine. The common disk file ''-REEDS, " Is then read, some 
local variables and display flags Initialized, and the exhaust cloud print header 
is written to the printing device. 

Useful Constant Terms for the Meteorological /Energy Factors . In 
preparation for the’ calculation of the stabilized height and cloud rise energy 
factor, two terms useful for the calculation at a meteorological factor for 
bouyant rise. 


GT 


9.8 


and 

.05129 (T - 273.16) 

QCj 

are calculated. 

Adiabatic or Stable Cloud Rise . Next, a loop is performed through 
the meteorological layers to simulate the cloud rise phenomenon, either 
adiabatic or stable. The criteria Is the argument of the Inverse cosine (ARCCOS) 
term in the bouyant cloud use equation. If the argument Is less than (*1), then 
the rise is "stable," and the flag so indicating, IAS, Is set to 1; otherwise, 
the rise Is "adiabatic," and IAS is 0. IAS Is Initialized in each layer to 1 
("stable" cloud rise). 
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Cloud Stabilization and Criteria for Surface Mixing Layer Height 
Selection . The slope (derivative with respect to altitude) of potential 
temperature (i.e., the **Lapse Rate"), the wind speed, and v/ind direction, in each 
layer are then calculated. The reader may recall the optional methods of 
calculating the lapse rate set up by the executive module, REEDS. If the least 
squares method (resulting in a lower cloud rise) is selected, a parameter IPTZ 
was set to 0. Otherwise, the alternate method of determining the lapse rate was 
chosen (IPTZ ■ 1), and is now calculated. In either case, a lapse rate is 
determined, and t''® energy factor calculated: 


a 




Atr 




where Z^ and Z^ are the current and initial layer altitudes, and Atj^ is the 
deposition time for the launch vehicle as defined by the vehicle's time versus 
altitude profile. The stability factor, STAB, is calculated from either lapse 
rate, determining adiabatic (STAB 10) or stable (STAB > 0, ll - H ctSTABi < -1) 
rise. If the 'rise is adiabatic, the rise time through the current layer is 
merely 

vT ; 


otherwise it is 
where * 


7T 


n -tan‘^(C3) 
/^STAB 




[1 - 1/2 g • STABl 
Vl - 11 - 1/2 a • STABl^ 


In addition to determining cloud rise, the three slopes calculated above become 
input factors to a deterministic calculation of the estimated surface mixing 
layer height to be performed after the conclusion of the rise loop and final 
stabilization height calculation: 

Stabilization Range and Transport Direction . As the cloud rises 
through the meteorological layers, the wind spef:d and direction change. This 
affects the cloud direction and resultant position, so these effects are 
determined. This results in a transport direction and a stabilization range 
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increment. At the conclusion of this loop, a stabilization range and height are 
calculated, later, the range and direction will be averaged over the rise to 
produce a net direction throughout the single surface mixing layer. 

Mixing Laver Height . Subroutine FTOP is then called to estimate the 
height of the surface mixing layer. To do this, the slopes of potential 
temperature, wind speed, and wind direction versus time are examined for 
"significant" changes. Just what is significant is somewhat subjective. 
Currently, three slope change magnitudes are stored in a data statement. They 
are: 


|i ^ 1 . 0.03 

I i ^ I ■ 0.015 





• 0.015 


for potential temperature, wind speed, and wind direction, respectively. When 
these are exceeded by the respective slope changes, the changes are considered 
"significant. " 

If no significant changes are found, the mixing layer height is set to 
twice the stabilization height. If only one of the variables has a significant 
change, the altitude of that change is selected. If there are more than one, 
they are sorted in descending order of magnitude. If there is an altitude at 
which all three variables have significant changes, this altitude is chosen. If, 
instead, there are two with coincident altitudes of significance, that altitude 
is chosen. If there are multiple multi-variable altitudes of significance, or 
several different altitudes which are non-coincident, then the largest change is 
selected. The following list summarizes the selection priority: 

1) The largest triple variable coincident significant 

change. 

2) The largest double variable coincident significant 

change. 
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3) The largest potential temperature significant change. 

4) The altitude of the single significant change. 

5) Twice the stabilization altitude. 

Wind Azimuth Uncertainty . A very Important driver of the turbulent 
diffusion process Is the wind azimuth uncertainty This variable Is 

computed by the subroutine “RS6AZ," u>^ng a theoretical and deterministic method 
developed by Panofsky from similarity theory. The alternative to a method such 
as this would be a statistical method. However, with soundings often used as 
Input meteorological sources, resolution at near-surface altitudes may not be 
good. This Is greatly Improved with tower data. Nevertheless, the similarity 
theory-based deterministic method Is an Important contribution for predictive 
application to single layer turbulent diffusion. 


The value of depends on the value selected from FTOP for the top of 
the surface mixing layer. That height is not final (i.e., it may be changed by 
the user if the Research or Operational mode is used, so final calculation of 
occurs when the mixing layer height is accepted. The form used for is 


T = - kf(B) 

Az - ln(z/zj - v(Rj) ’ 


CTV 


U 

k 

f(8) 


z 


standard deviation of lateral component of turbulence 
(m/sec) 

mean wind speed (m/sec) 

Von Karman's constant (».4) 

Piecewise linear fit to experimentally measured values 
of av/^ as a function of the non-dimensional stability 
ratio, B 

altitude (m/ 


Zq =* roughness length (m) 

■^(RI) = a universal function of the Richardson number. 


2-26 


J:/ 


This relation was derived by Lettau (unpublished) and Panofsky (1963), serving 
both stable and adiabatic atmospheres (International Journal of Air and Water 
Pollution, Vol. 9, 1965, Panofsky and Prusad on ''Similarity Theory and 
Diffusion"). The non-dimensional stability ratio, B, Is given by 


B - 



A(Pj) 

^z 


where 


9 

7 

u 

T 

dP/ 


■Iz 


2 

gravitational acceleration (m/sec ) 

the geometric mean of layer top and bottom (m) 

mean wind speed (this layer) (m/sec) 

absolute temperature (°K) 

lapse rate (vertical gradient of potential 
temperature). 


Depending on the value of B, f(B) Is selected from 
f(B) ■ a + b (B + c), where 


1 A*2.7, b»c*o for B < -.008 

» 2.7, b - 112, c ■ .008 -.0081 B < -.00175 

- 3.4, b * 725.5, c » .00175 -.00175 < B< .0008 

» 1.55, b * 38.04, c « -.0008 .0008 < B < .029 

- 2.35, b » 5.43, c - .029 .029 1 B 

The Richardson number for 'p(RI) cannot be calculated from its definition 
(normalized lapse rate per unit vertical wind speed gradient squared), because 
sufficiently accurate measurements of the vertical wind speed gradient do not 
normally exist. It can be transcendental ly obtained from: 


R 


I 


pn(z/zQ) - -^(Rj)"!^ 

L J 


where (RI) « (1-16 Rj^)”^, unstable conditions 

(1-7RI)"^, stable conditions 
(another universal function). 
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and 


2 

’^/(RI) ■ In + In -^^1“ + j - 2 Tan”^(x), unstable condition 
■ 7 Rj(l-7RI)"^. stable conditions 


where X - (1-16RI)^. 

For stable conditions, 

sjTi - (-1 + JTTTT^ )/2Aj 

where Aj^ * 7 Bk, and A^ ■ - B (k-1), for k ■ In l/l^- 

For unstable conditions, 

— 5 r-,.B . 0 

16 XMln z/Zq + .50864 - 2 In (1+x) - ln(l+x^) + 2tan'^xl ■ 

is solved transcendental ly for X, which yields RI. 

In RSGAZ, the transcendental function statement and its necessary 
derivative precede the active statements (follow the detailed flow chart in 
Figure 2-5). The critical data (transferred via common) is selected from the 
input layers as a result of the input arguments passed. The potential 
temperature gradient (Lapse Rate) is then calculated at the indicated altitudes, 
followed by the Richardson number for either stable or unstable conditions. The 
appropriate value of f(B) is then obtained for the final calculation of 
value of returned is constrained between the limits of 7 degrees and 30 
degrees, to keep resultant uncertainties physically realizable. 

Interactive Alteration of the Surface Mixing Layer Height and the Wind 
Azimuth Uncertainty . Once the mixing layer height and subsequent wind azimuth 
uncertainty are estimated by the program, the user is given the opportunity to 
change each, if the RESEARCH or OPERATIONAL mode has been chosen. In the 

production mode, the user must enter both. The user is first prompted to use the 
estimated layer height or an altered value. If the altered^^iajue choice is 
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Figure 2-5. Detailed Structure Flow Diagram of Subroutine RSGAZ 







ORIGINAL PAGE IS 
OF POOR QUALITY 




Figure 2-5. Detailed Structure Flow Diagram of Subroutine RSGAZ (continued) 
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Figure 2-5. Detailed Structure Flow Diagram of Subroutine RSGAZ (continued) 
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Figure 2-5. Detailed Structure Flow Oaigrain of Subroutine RS6AZ (continued) 







Figure 2-5. Detailed Structure Flow Diagram of Subroutine RSGAZ (continued) 
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Figure 2-5. Detailed Structure Flow Diagram of Subroutine RSGAZ (continued) 
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Figure 2-5. Detailed Structure Flow Diagram of Subroutine RSGAZ (continued) 
















selected (via entry of an Integer > 8) an interactive selection of the layer 
height Is assumed. In either case, the coimon data are written to **"REE0S," and 
the meteorological data plotting routine "RMTPS" Is called. Upon conclusion of 
RMTPS, control Is returned to RCcDS. If met plots only were desired, the call to 
'•RMTPS'* Is followed by a program STOP for Immediate termination of execution. 
Otherwise, the height selected during "RMTPS" becomes the estimated mixing layer 
height for the balance of the calculations. 

Contaminant Source Strength . Calculation of the mixing layer height 
determines the source strength of HCL (from which the other contaminant strengths 
are determined) In that It determines the maximum altitude over which the rocket 
vehicle effluent production algorithm may operate; contaminants dumped above the 
surface mixing layer are not permitted to diffuse Into the layer. 


Final Transport Speeds Direction, Wind Azimuth Uncertainly and 
Stabilization. . Definition of the mixing layer height also permits correct 
averaging of transport speed and direction over the layers within the mixing 
layer, so this Is accomplished and the results written on the selected printing 
device. The user Is then given the chance to change the wind azimuth 
uncertainty. If the RESEARCH or OPERATIONAL mode is active. If not "N" for No, 
the user enters the uncertainty desired in wind azimuth (degrees) and is 
printed on the printing device. Since the mixing layer height selected (Hj) may 
have resulted in part of the cloud being above the mixing layer top, and 
therefore being discarded, the "effective" stabilization heat (H^g) of a 
"revised" cloud containing only that portion within the mixing layer must be 
calculated. This is done according to 


H 


SE 



Hj - < R|. 

I«T - ''c 


where is the stabilization height and R^, the cloud radius. 
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I 



where ^Xo ^ s the initial cloud radius. 

The stabilization height is written to the printing device, the coimion 
data are written to the disk file, and the effluent concentration module ''RCONS*' 
is called. Following this, a program STOP is generated, returning control to the 
"REEDS" executive module. 

2.1.3 The RMTPS Meteorological Data Plotting Module 

.The "RMTPS" module generates a plot of the Meteorological Profile 
Variables of wind speed and wind direction as a function of altitude, and the 
initial and effective stabilized cloud shape. The mixing layer is also indicated 
on the olot and can be interactively determined from user input. As an 
additioni' feature, the sounding (or other) sampling points are Indicated as 
additional tic marks on the ordinate altitude scale. The "RMTPS" module also 
will draw the met plot forms if the user has selected that option in "REEDS." 

Forms Generation . After the common disk file is read to load the 
common data, the RMTPS module checks for the forms switch, MFORM. If on (MFORM * 
1), the form is generated; the first two hundred active lines (or so) are 
dedicated to this purpose. 

Met Plots . After the forms are generated, or when forms generation 
is bypassed, the MONLY common parameter (OPTION control pseudo-block) is checked 
to see if MET plots are skipped altogether (or in a forms generation run or in a 
repeat run designed to prevent needless replotting of the MET-related data). The 
wind direction data Is rotated to conform to the preset site orientation, and is 
scaled to plot units. The launch data and time, site data, surface pressure and 
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density, stabilization height and wind azimuth uncertainty are all encoded into 
ASCII character strings and written by the plotter onto the MET plot. This is 
followed by the dry temperature, potential temperature, wind speed, and wind 
direction (all surface values), which are likewise encoded and written by the 
plotter. In each case, the encoding is accomplished by the combination of "CODE" 
calls and subsequent WRITE to the character string, lALPHA. The character string 
is then written at the designated plotter location via the "CHAR" subroutine. 

Curve Drawing . Actual curve drawing then proceeds by pulling the 
array values for the various arrays out and using them as data for the "CHAR" 
calls. This results in the indicated symbol from the symbol array ICRVT being 
plotted. In the case of the line plots, either "LINE" (for solid lines) or 
"OLINE" for dashed lines are used for curve generation on the plotter. 

Ax is /Tic Marks. Curve drawing is then followed by axis generation. 
This is accomplished with simpler calls to LINE, whether for axes or for tic 
marks. 

Cloud Drawing . The cloud schematic is drawn by the "CLOUD" 
subroutine in the lower right quadrant of the plot. Relative distances with 
respect to mixing layer are kept to approximate scales. 

Mixing Layer Height . At the appropriate height on the altitude scale 
(ordinate), the estimated mixing layer is drawn with the "LINE" routine. The 
"MOVEM" subroutine is then called to prompt the user to accept or replace the 
indicated height. The choices are; 

"CO" for continue 

"UP" for move up one met level 

"DO" for move down one met level 

After each prompting, the new line is drawn in. When it is finally accepted by 
entering "CO," the line is labeled by the encoded layer height (meters). 

Effective Cloud . Some cloud parameters (radius, height, etc.), 
calculated in RCLOS but not saved, are recalculated. This permits redrawing the 
cloud schematic with the "effective" parameter. This time, however, the "LINE" 
routine is used with the cloud extent data to produce a shaded ellipse 
superimposed over the original cloud outline. After this operation, the common 
data are written to the disk file and a program "STOP" is generated, returning 
control to "RCLDS." 
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2.1.4 


The RCONS Centerline Concentration Module 


The RCONS concentration module computes the range and azimuth of the 
cloud center relative to the pad, calculates and plots the centerline dose and 
concentration of the major hazardous effluents, computes the off-centerline dose 
and concentration for specific locations, and calls the "RISOS" dosage and 
concentration isopleth plotting module. The ''RCONS'V'RISOS" combination 
comprises the "kinematic phase" of the effluent diffusion problem. Refer to 
Figure 2-3 for the RCONS top-level flow. 

The turbulent diffusion problem may be described by the following 
diffusion equation: 

2t * ^x(r.t) * V • [ K(r,t,p,T) • Vx(r,t)] 

where (r,t) is the scaler concentration of the diffusing effluents, <V> the 
expected value of the wind velocity and K (r,t,P,T) the diffusion diagonal 
tensor which is a function of temperature, pressure, position and time. The 
separation of the problem into two (2) distinct phases, a thermodynamic phase 
resulting in approximate equilibrium with the ambient environment and a 
kinematic phase, with diffusion from the thermodynamically stable state, 
facilitates the linearization of the equation and subsequent separation of 
variables solution. We have already looked at the first phase (RCLDS), and now 
wish to examine the latter. 

The first phase left the exhaust cloud initially spherical in shape 
with the cloud center at a stabilized altitude. Then, the revision of the cloud 
occurred to account for surface mixing layer effects. Now, the ellipsoidal 
expansion of the revised cloud occurs. The generalized dosage equation driving 
this expansion is composed of four terms: 

nncACF . f 'I « /lateral) ^ /vERTICAl) ^ /dEPLETIOn) 

DOSAGE V DOSAGE/ * V TERM / \ TERM / ^ \ TERM / 

In subsequent modeling, the depletion of material from the cloud by gravitational 
settling and rain scavenging is included, but in the current REEDS model, no 
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depletion of this type 1s Included. Therefore, the driving dosage algorithm Is 


0 |x, 

+ exp 

+ exp 
+ exp 


y. Zg < z < Zt( 
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■- (H -2Zg + z)*' 
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2a, 


r-[2f (Zj . Zj) t {H - z)]=l r-[2i (^T - ^8^ 

L 257 J * L- 5 S 7 ~ 

|- -[2i (z^ - z^)*^ (h - 2Ze ^ z)]» j 


(H - 



where : 

Qm » Source Strength (total mass In the indicated surface 
layer) 

T,B subscripts indicate top and bottom of stabilized cloud, 
respectively. 

H » Height of Cloud Centroid 

Z » Height within layer (Zgl z 1 z^) 

cr , <j * Standard Deviation of the Time Dependent Cloud 

^ Dimensions in the y and z directions, respectively. 

Centerline values are obtained at y - 0, while surface values are obtained at z - 
0. From the dosage, the concentration is formed by the expanding front as 




P(Xi'y>z,t) u 
'^x 
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where U is the average wind speed, and is the along-wind cloud dimension. The 
balance of the kinematic phase modules are concerned with establishing positions 
and timings relative to the cloud center and time of maximum 
dosage/concentration, and the relative distribution of dosage/concentrations 
among the effluent constituents. 

Centerline Dosages and Concentrations . After some input terminal 
initialization, the standard coninon disk file read, and some useful parameter 
initialization, the range down the cloud centerline is computed and expressed 
relative to the pad. The dosage algorithm "guts" exists in subroutine OFEXP, 
where the various sums in the above equation are computed. DFEXP is called with y 
■ 0 off centerline distance to obtain the centerline maximum dose and 
concentration. The arguments of DFEXP are the surface mixing layer height and 
the concentration desired. DFEXP can then obtain the appropriate dimensions 
corresponding to that concentration. However, the specification of the number 
1000 for concentration is a key that causes DFEXP to produce only values at the 
indicated range, and report the standard deviation of the cloud dimensions in the 
transverse direction. 

Cloud Lenq%i and Arrival /Departure Time . Based on the calculations 
of cloud extent made by DFEXP, the arrival and departure times at the indicated 
range are computed and printed with the range and azimuth relative to the pad, at 
250 meter increments along the cloud centerline. 

RCONS User Prompts . Most prompts for RCONS user action are stored in 
the RCONS question file "7RC0NS," and are accessed via the DREAD subroutine. 

Centerline Plots . If the current operational mode is not production, 
the user is then prompted to decide whether or not to produce actual plots of the 
centerline dosage and concentration. If the response is not "N" for no, 
subroutine CPLOT.is called with the centerline range array generated above. The 
plot is generated with a pre-set scale factor in order to coincide with 
traditionally generated forms. The default response to the centerline plot 
prompt is "YES." 

Maximum Dosa^a/Concentration . The maximum dosage and concentration 
will occur along the cloud centerline. The point at ground level where this 
occurs had been calculated and saved during the above centerline dosage/concen 
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tratlon calculation. They are now written out on the printing device selected 
with pad-relative range and azimuth. 

Off-Centerline Concentrations /Dosages . If the current operating 

mode is PRODUCTION, the run is effectively complete and transfer of control is to 
the end of RCONS, and ultimately to the REEDS executive module for termination or 
restart. If the current mode is not PRODUCTIDN, however, the user is prompted 
for "off -center line dosages/concentrations. “ If the answer is not "N" for "No," 
off -center line concentrations/dosages are assumed to be desired (which is the 
default case), and the ORIGIN subroutine is called to set the origin for the 
desired launch pad and map. Application of the REEDS code to other sites can be 
accomplished by appropriately modifying the vehicle/pad association data in the 
ORIGIN subroutine. From ORIGIN, the user is prompted for the pad number and 
land/sea map combination. The choices are listed by the prompt, depending on the 
vehicle selected during REEDS input (certain vehicles imply certain pads). The 
user then inputs the code listed in the prompt which corresponds to the pad/map 
combination he desires. The user will now be prompted for pad-relative off- 
centerline ranges and azimuth. However, if the operating mode is RESEARCH, the 
user is first prompted (from the programl to enter a cloud centerline 
(stabilization-relative) range. This is then transformed to a pad-relative 
range and azimuth using the setup from ORIGIN. This is an aid in the 

verification of data from the centerline calculations and/or an additional 
centerline/off -centerl ine consistency check. In either RESEARCH or OPERATIONAL 
modes, the user is then prompted for the pad-relative range followed by the pad- 
relative azimuth. The entered range and azimuth are printed on the terminal 
display, and then on the printing device. The coordinates are plotted on the 
Isopleth plotting map requested by ORIGIN according to the setup resulting from 
the requested pad/map combination. The point is labeled by the ordinate number 
of the range/azimuth pair .(first entry » 1, second entry » 2, etc.). Two 

additional points are then computed at the requested range, but the first at 10° 

less azimuth and the second at 10° more azimuth. All three range/azimuth pairs 
are then transformed to cloud center -relative coordinates, and DFEXP is called. 
The concentration and dosage is evaluated at each point and the +10° data used to 
produce an uncertainty in the concentration/dosage for each effluent. The 

resulting data is printed on the printing device. The format is: 
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MATERIAL 


CONCENTRATION ± CONCENTRATION DOSAGE + DOSAGE 
(PPM) UNCERTAINTY (PPM-SEC) UNCERTAINTY 

(PPM-SEC) 

The above off-center! ine calculation is repeated until a zero (0) range is 
entered. 

Isopleth Plotting . The user is then prompted for Isopleth plotting. 
If the answer is not "N“ for "No," the Isopleth plotting module is scheduled, 
otherwise RCONS is terminated and control returns to REEDS through RCLOS, for 
program termination/restart. 

2.1.5 The RISOS Isopleth Plotting Module 

The RISOS module essentially repeats off -center 1 ine calculations using 
OFEXP, except it uses a clever series of calls resulting in iso-concentration 
contours. Refer to Figure 2.3 for RISOS top level flow. 

Isopleth Initialization . The user terminal is first initialized and 
the common data read from disk via subroutine HWDIS. The ORIGIN subroutine is 
called again in case the call in RCONS was bypassed by rejection of the off- 
centerline option. This makes sure the plotting origin is established relative 
to the vehicle pad/map configuration desired. * 

Cloud Rise Ground Track . The cloud rise ground track from the pad to 
stabilization is then plotted by calculating the range of the cloud centroid from 
the time required to rise through each meteorological layer. The times for each 
layer were stored during the cloud rise module (RCLDS) execution. The average 
wind velocity in each layer is integrated over the stored time to provide the 
incremental range that is plotted. Since the wind velocity is measured "from" 
the source, the cloud transport direction must use a factor of ^ to properly 
orient its movement. The process is terminated when the layer containing the 
stabilized cloud center is reached. However, an additional segment is required, 
because of the use of the "effective" stabilization height, which replaced the 
stabilization height when placement of the surface mixing layer height severed 
part of the cloud (for diffusion purposes). The stabilization point is then 
labeled with a "+." Next, the point of maximum concentration at ground level, 
stored during RCONS, is labeled with an Finally, the cloud track downwind 
from stabilization is plotted with a straight line using the transport direction 
from that point. 
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3.0 REEDS OPERATION 


The basic form of storage for REEDS is in the form of source code files 
for the program/subroutine packages and ASCII data files for the input data. 
Ultimately, the source code files are compiled and thereby translated into 
machine language modules. The latter are loaded, SP'd, OF'd and RP'd (see 
terminal users manual). The ''RU“ command is then used to execute the RP'd 
modules while arguments are used to pass special information to the program. The 
compiled, loaded, and SP'd files are convenient to use, since they merely have to 
be RP'd prior to execution— a very quick process. One drawback is that there is 
little way of verifying that the RP'd module is the result of the source file as 
it appears. Also, since source files must be saved anyway, the existence of a 
large compiled file may seem redundant, especially if not used often. 
Nevertheless, because of the lengthy loading process, the compiled files are 
often stored. Thus, for purposes of showing a run, the compiled files are 
assumed to exist in the RP'd state . Maintenance files have been built which 
automatically perform the operations leading to an RP'd program. It is the 
purpose of this section to explain actual REEDS execution. 

3.1 THE TEST CASE 

A test case was chosen, using the space shuttle as the vehicle and the 
8 A.M. EDT 16 November 1969 rawinsonde as the meteorology. The file name of the 
latter data was changed to DATAOl, just for convenience, and can be seen in 
Section 4.4 of this report. The CRT terminal used was equipped with a hardcopy 
device, making it possible to get intermittent copies of the screen to aid the 
narration. 

Run Time Parameters . The REEDS system was designed to serve MSFC 
research needs and KSC operational needs. One difference between these is that 
KSC has three plotters, each one assigned to a major plot output. MSFC has one 
plotter for all three plots. The allocation in the run statement of printing and 
plotting devices to logical units makes this possible. The user enters: 
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Isopleth plotting 
Device (7922 
Plotter selected) 


RU, REEDS, 1, 6, 


12 , 12 , 12 



Input Device 
(Users CRT-See 
Terminal Users 
Guide) 


y 



Printing Device 
(Line printer 
selected) 


Centerline Plotting 
Device ^ 
Meteorological (7922 plotter selected) 
Plotting Device 
(7922 plotter 
selected) 



At this time the screen appears as In Figure 3-1. The default responses are 
Indicated on the CRT screen when they appear by easily noted enhanced characters. 

Operational Configuration . The next five (5) entries deal with 
configuring the logical path through the modules. The first prompt asks about 
the low altitude curve fit. We choose not to use the least squares fit for the 
potential temperature gradient and answer "No." We have renamed the 16 Nov 1969 
rawinsonde file DATAOl for this run, and therefore select the default, "YES," 
response for the meteorological data file name (response-2). For the third 
Input, the default^ OPERATIONAL mode Is chosen, followed by a similar choice for 
Met Plot Form rejection (4th) and Met plot creation (5th). ' The printer output 
default (“Yes") is selected, and the screen appears as in Figure 3-2. 

Launch Specific Parameters . The default launch time is the current 
time on the computer's system clock. It was decided to select the time rather 
than enter either a realistic or other time. Likewise, the default vehicle 
(space shuttle) is chosen as the launch vehicle. At this point, the screen 
appears as in Figure 3-3. 

Initial Printed Output . At this point, the initial printer data 
becomes available, printing the header and meteorological sounding data as shown 
in Figure 3-4. The printout reflects all input entered to date, except the Met. 
Plot Form/Met. Plot configurational data. 

Estimated Mixing Layer Height . The RCLOS cloud rise module has 
estimated the height of the top of the surface mixing layer, printing 2552 meters 
on the screen. The user is then prompted to raise or lower this height. The 
cloud rise has been calculated and stored up to this height, and RMTPS has been 
called. 
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RU, REEDS, 1,6, 12, 12, 12 

*****NASA/MSFC EXPERIMENTAL MULTILAYER 0IFFU15 AUG 1981**** 


Figure 3-1. Code Initiation 


*****NASA/MSFC EXPERIMENTAL MULTILAYER 0IFFU15 AUG 1981**** 
DO YOU WANT THE LOW ALT. CURVE FIT? YES NO NO 


LOW ALTITUDE CURVE FIT NO 

USE DEFAULT COMMON DATA FILE (DATA01) ? YES 

RUN TYPE: OPERATIONAL 

GENERATE THE MET PLOT FORM? NO 

MET PLOTS? YES 

GENERATE PRINTER OUTPUT? YES 

DEFAULT LAUNCH TIME/DATE? 1833 CDT 15 AUG 1981 YES NO 


Figure 3-2. Configurational Paths Added 
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★★★★★NASA/MSFC experimental multilayer 0IFFU15 AUG 1981*^** 


00 YOU WANT THE LOW ALT. CURVE FIT? YES NO NO 

LOW ALTITUOE CURVE FIT NO 

USE OEFAULT COMMON OATA FILE (OATA01) ? YES 

RUN nPE: OPERATIONAL 

GENERATE THE MET PLOT FORM? NO 

MET PLOTS? YES 

GENERATE PRINTER OUTPUT? YES 

OEFAULT LAUNCH TIME/OATE? 1833 COT 15 AUG 1981 YES NO 

SELECT LAUNCH VEHICLE: • ^ SHUTTLE-S 

ESTIMATED TOP OF SURFACE LAYER(M) IS: 2552.4 


MOVE TOP OF SURFACE LAYER: UP DOWN CONTINUE 


Figure 3-3. Launch Specific Data Added 
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## ######### ifr ### A*# #'*#### ir 

isifititifitiei^’kirii kirirititirit^iikk 

#★★★##★★*★★ NASA/MSFC REED CODE USING SINGLE LAYER MODEL 21 DEC 1979 ★*★★★★★★★★★ 

kki^itififirif^ir #########★★ 

★★★★★★★★★★★ ★★★★★★★*★★★ 

★★★★★★★★★★★ REEDS irkiritkkiritititk 

’kiirk’kititititifitit ★★★★★★★★★★★ 

kkirkiritifirkkkkkifififirkiriritkirkitiriritiritiririfiririfkirifkitiririritiritiririfitiriririritirir^kiririrififitirkifiritifitkiritiririrititifir 


RUN 1 USING DATA FILE DATAOl 
SHUTTLE LAUNCH VEHICLE 

PREDICTION TIME; 1007 CDT DATE: 24 AUG 1981 

DATA FILE LEADER INFORMATION: 

TEST NBR 03717 04834 0-24HR 
RAWINSONDE RUN AN/GMO-4 
CAPE CANAVERAL AFS, FLORIDA 
ASCENT NBR 0543 

TIME: 800 EDT DATE: 16 NOV 1969 


SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
SSSSSSSSSSSSSSSSSSSSS SSSSSSSSSSSSSSSSSSSSS 
SSSSSSSSSSSSSSSSSSSSS SOUNDING SSSSSSSSSSSSSSSSSSSSS 
SSSSSSSSSSSSSSSSSSSSS SSSSSSSSSSSSSSSSSSSSS 
SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS5 


SURFACE DENSITY (GM/M**3): 1273.00 


LAYER 

ALTITUDE 

DIRECTION 

SPEED 

TEMP 

POT-TEMP 

0 P TEMP 

PRESSURE 

NO. 

(FEET) 

(METERS) 

(DEGREES) 

(M/SEC) 

(DEGREES CENTRIGRADE) 

(MILLIBARS) 

1 

16 

5 

220 

.4 

7.1 

4.96 

23.2 

1027. 

2 

735 

224 

10 

.8 

9.6 

9.60 

22.9 

1000. 

3 

754 

230 

10 

.8 

0.7 

9.76 

21.7 

999.3 

4 

2126 

648 

27 

.7 

6.5 

10.66 

16.6 

950.0 

5 

2172 

662 

28 

.7 

6.4 

10.69 

15.3 

948.5 

6 

3579 

1091 

7 

.4 

6.3 

14.91 

13.6 

900.0 

7 

3629 

1106 

6 

.3 

6.3 

15.04 

13.6 

898.6 

8 

5121 

1561 

316 

.3 

7.2 

20.63 

10.6 

850.0 

9 

5960 

1817 

335 

.2 

7.7 

23.78 

8.7 

824.2 

10 

6754 

2059 

9 

.2 

6.4 

24.96 

10.7 

800.0 

11 

8374 

2552 

41 

.1 

4.9 

28.51 

7,6 

753.5 

12 

3485 

2586 

41 

.1 

4.9 

28.92 

3.5 

750.0 


Figure 3-4. Initial Printer Data-Header and Meteorological Data 
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Initial Plotter Output . The meteorological data are plotted, since 
^«)t. Plots were selected; however, the form Is omitted because a blank form 
already exists and Is used. The blank form appears In Figure 3-5. After the 
Meteorological data are p jtted, the plot appears as In Figure 3-6. The 
estimated mixing layer top Is drawn, but not labeled, and the meteorological 
sounding data points are Indicated by the additional Irregularly spaced tic marks 
on the left ordinate. The program Is left In an Input mode, awaiting 
Instructions as to mixing layer height modification. 

Mixing Layer Height Modification . As a demonstration of the program 
options rather than a reflection of reality, the cloud Is shown stabilized at 871 
meters, and the mixing layer height Is dropped four (4) times to 1106 meters. As 
can be seen In Figure 3-7, this Is accomplished bj# entering four (4) consecutive 
"D0'“s for (DOWN) and finally a "CO" (fc** continue). 

Final Cloud Rise . With the mixing layer finalized, the cloud rise Is 
modified to reflect "severing" of part of the cloud by the boundary and 
subsequent recalculation of the stabilization height (to become the "effective" 
height) and total rise time. The resulting "effective" cloud Is shown 1n Figure 
3-8, and the cloud rise summary data printed on the printer Is shown In Figure 3- 
9. 

Final Wind Azimuth Uncertainty . The wind azimuth uncertainty is 
calculated from the mixirki layer height and printed on the terminal display 
(Figure 3-7). The user Is given a chance to modify this value. As the figure 
shows, the estimated value was accepted. 

Centerline Concentrations and Dosages . The table of calculated 
centerline HCL dosages and concentrations are printed as shown In Figure 3-10. 
The user is prompted for centerline plots of concentration and dosages. These 
are desired (Figure 3-11) so a blank centerline form is loaded onto the plotter 
(Figure 3-12). The centerline dosage and concentration are plotted on the 
plotter, and the plot is labeled (Figure 3-13). The maximum ground concentration 
and dosage are then printed on the printer, as are the pad-relative coordinate 
(Figure 3-14). All dosages and concentrations are HCl dosages and concentra- 
tions. All other effluents are determined from HCl, the chief effluent. 
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Figure 3-5. Blank Meteorological Form 
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Figure 3-6. Meteorological and Initial Cloud Data Plotted with 
Estimated Mixing Layer Height 
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USE A SIGMA FOR WIND AZIMUTH ANGLE OF; 

IS: 

30,0 

1106.12 

YES 

NO 

USE A SIGMA FOR WIND AZIMUTH ANGLE OF: 

30.0 

YES 


CENTERLTNE concentration PLOT DESIRED ? 


YES 

NO 


Figure 3-7. Modifying the Mixing Layer Height 
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Figure 3-8. Modified Cloud - Showing Mixing Layer Height Effects 


LEVEL 


EXHAUST CLOUD 


ALTITUDE 


RISE TIME 

RANGE 

DIRECTION 

(METERS) 


(SECONDS) 

(METERS) 

(DEGREES) 

224.0 

STABLE 

7.6 

4.7 

165.0 

229.8 

STABLE 

8.0 

5.0 

166.5 

648.0 

STABLE 

53.3 

39.3 

194.6 

662.0 

STABLE 

55.5 

40.8 

195.1 


★★★★CLOUD STABILIZATION**** 

HEIGHT(M): 1004.2 

STABILIZATION TIME AFTER LAUNCH(SEC): 174.0 
RANGE FROM PAD(M): 111.7 

DIRECTION FROM PAD(DEG): 195.3 

PRESSURES ARE: 1027.00 1000.00 999.30 

TBHT « .733812E+03 


****TOP OF SURFACE LAYER METEOROLOGICAL PARAMETERS**** 
HEIGHT(M):- 1106.1 
WIND DIRECTION(DEG): 6. 

WIND SPEED(M/SEC): .3 


****DIFFUSION PARAMETERS**** 

MEAN SPEED(M/SEC): .7 

MEAN TRANSPORT DIRECTION(DEG) : 191.3 


SIGMA OF WIND AZIMUTH ANGLE, SIGAZ: 30.0 


EFFECTIVE CLOUD HEIGHT (M): 733.8 

Figure 3-9. Final Exhaust Cloud Rise Data 
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Figure 3-10. Centerline Hcl Concentratlons/Oosages with Cloud Arrival /Departure 
Times and Pad-Relative Positions 


DO 


DO 

DO 

DO 

CO 


CALCULATED TOP OF SURFACE TRANSPORT LAYER 
USE A SIGMA FOR WIND AZIMUTH ANGLE OF: 

IS: 

30.0 

1106.12 

YES 

NO 

USE A SIGMA FOR WIND AZIMUTH ANGLE OF: 

30.0 

YES 


CENTERLINE CONCENTRATION. PLOT DESIRED ? 


YES 

NO 

YES 

CENTERLINE CONCENTRATION PLOT DESIRED ? 


YES 



Figure 3-11. Centerline Concentration/Dosage Plots Added 
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Figure 3-12. Blank Form for Centerline Concentration/Dosage Plots 



Figure 3-13. Centerline Dosage/ Concentration Plots 
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Off -Center line Concentration/Dosaqes . The user is thon prompted for 
off -center line concentrations/dosages. These are selected, causing a call to 
the ORIGIN subroutine and a further prompt for the desired pad number/map type. 
Of the choices listed, pad 39 was chosen, and a sea map (Figure 3-15) rather than 
a land map (Figure 3-16 is chosen. As will be seen, this is a mistake, as better 
resolution and coverage would result for a land map. 

The user is then prompted for range and azimuth (pad-relative) to the 
point for which off -center line calculations are desired. The prompting 
continues (Figure 3-17) until a zero (0) range is entered. During this time, the 
off-center! ine effluent concentrations, their associated +10° (in azimuth) 
uncertainties, and dosages and their respective uncertainties, are printed on 
the printer (Figure 3-18). The points are also plotted on the plotter, and 
labeled with their ordinate index (1st one « 1, 2nd * 2, etc.) (Figure 3.19). 

Isopleth Plotting . The user is then prompted for isopleth plotting. 
This is accepted and the default set of concentration levels is listed on the 
terminal display, and the user is asked if he intends to use these, or input his 
own (Figure 3-20). These are accepted, and the isopleths are drawn and (Figure 
3-21) labeled. 

Clouu Arrival /Departure Time . Finally, the cloud position and extent 
as a function of time are printed in a concluding table (Figure 3-22). 

Termination . The user is then asked for termination ("T") or restart 
("R") termination is chosen (Figure 3-20). 
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CONCCNTHA nON VAI U£S 
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COMCf NTRAT!OM VALUES Wm) 










SPECIFIC SITE PREDICTIONS (PLOTTED)? 


YES 


NO 


39. SEA MAP 

ENTER RANGE(M) AND AZIMUTH(DEG) : (0 TERMINATES PROCEDURE) 

RANGE: 1000 AZIMUTH: 192 

RANGE: AZIMUTH 


ENTER RANGE(M) AND AZIMUTH(DEG) : (0 TERMINATES PROCEDURE) 


RANGE 

1000 

AZIMUTH 

192 

RANGE 

2500 

AZIMUTH 

192 

RANGE 

5000 

AZIMUTH 

192 

RANGE 

7000 

AZIMUTH 

192 

RANGE 

mo 

AZIMUTH 

192 

RANGE 

10000 

AZIMUTH 

192 

RANGE 

12000 

AZIMUTH 

192 

RANGE 

15000 

AZIMUTH 

192 

RANGE 

20000 

AZIMUTH 

192 

RANGE 

25000 

AZIMUTH 

192 

RANGE 

0 

AZIMUTH 


RANGE : 

: 0 

AZIMUTH: 


CONCENTRAflON ISOPLETH PLOT DESIRED? 


YES NO 


Figure 3-17. Specific Site Range/Azimuth Entry 
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♦♦♦♦CONCENTRATIONS AND DOSAGES WITH 10 DEGREE UNCERTAINTIES^**^ 


RANGE(M): 5000.0 




AZIMUTH (DEG): 

170.0 




MATERIAL 

CONCENTRATION(PPM) 


DOSAGE 

(PPM) 

HCL 

.184E+01 +/- 

.050 

’,09.108 +/- 

57.684 

CO 

.676E-02 +/- 

.000 

7.730 +/- 

.211 

C02 

.403E+01 +/- 

.110 

4606.357 •►/- 

125.984 

AL203 

.473E+01 +/- 

.129 

5409.617 V- 

147.953 

NO 

.322E-02 +/- 

.000 

3.681 +/- 

.101 

RANGE (M): 5000.0 




AZIMUTH (DEG): 

190.0 




MATERIAL 

CONCENTRATION(PPM) 


DOSAGE 

(PPM) 

HCL 

.306E+01 +/- 

.285 

3494.906 +/“ 

326.191 

CO 

.112E-01 +/- 

.001 

12,809 +/- 

1.195 

C02 

.668E*01 +/- 

.623 

7632.983 +/- 

712.411 

AL203 

.784E-HD1 +/- 

.732 

8964.029 V- 

836.641 

NO 

.534E-02 +/- 

.000 

6.099 +/- 

.569 


RANGE(M): 7500.0 




AZIMUTH (DEG): 

190.0 



(PPM) 

MATERIAL 

CONCENTRATION(PPM) 


DOSAGE 

HCL 

.212E+01 +/- 

.219 

2428.497 +/- 

250.218 

CO 

.779E-02 +/» 

.001 

8.900 +/- 

.917 

C02 

.464E+01 +/- 

.478 

5303.912 +/- 

546.484 

AL203 

.545E+01 +/- 

.561 

6228.813 +/- 

641.781 

NO 

.371E-02 +/- 

.000 

4.238 +/- 

.437 

RANGE(M): 10000.0 




AZIMUTH (DEG): 

190.0 




MATERIAL 

CONCENTRATION(PPM) 


DOSAGE 

(PPM) 

HCL 

.163E+01 +/- 

.176 

1859.782 +/- 

201.469 

CO 

.596E-02 +/- 

.001 

6.816 +/- 

.738 

C02 

.355E+01 +/- 

.385 

4061.821 +/- 

440.015 

AL203 

.417E+01 +/- 

.452 

4770.125 +/- 

516.745 

NO 

.284E-02 +/-+/- 

.000 

3.246+/- 

.352 

RANGE(M): 12500.0 




AZIMUTH (DEG): 

190.0 




MATERIAL 

CONCENTRATION(PPM) 


OOSAGE 

(PPM) 

HCL 

.132E+01 +/- 

.147 

1506.710 +/- 

168.214 

CO 

.483E-02 +/- 

.001 

5.522 +/- 

.616 

C02 

.288E+01 +/“ 

.321 

3290.701 +/- 

367.386 

AL203 

.338E+01 +/- 

.377 

3864.536 +/- 

431.451 

NO 

.230E-02 +/-+/- 

.000 

2.630 +/- 

.294 

Figure 

3-18. Off-Centerl ine 

Concentrations and Dosaaes 

with +10° 

Uncertainties for aV 

1 Effluents 
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CONCENTRATION VALUES (ppm) 



Figure 3-19. Off-Centerline Point Plotting 








CONCENTRATION ISOPLETH PLOT OESIRED? YES NO 

YES 

DEFAULT ISOPLETH CONCENTRATION VALUES ARE: CONC(l) • .515 

CONC(2)-2.573 C0NC(3)»3.859 CONC(4)»1.0 

USE DEFAULT ISOPLETH CONCENTRATION VALUES ? YES NO 

YES 

T 

***** PROGRAM REEDS HAS TERMINATED NORMALLY **** 

REE05 : STOP 0000 

Figure 3-20. Default Isopleth Values, Isopleth Plotting, 
and Program Termination 


» 
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Figure 3-21. Isopleth Plotting 


CLOUD LOCATION AND DIMENSIONS 

TIME FROM LAUNCH RANGE AZIMUTH DIAMETERS (METERS) 
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Figure 3-22. Cloud Position and Extent as a Function of Time 


SECTION 4.0 LISTINGS 


This section contains the current listings of critical elements of the 
REEDS system. They are included as follows: 

SECTION 4.1 Main Program Modules and Subroutines Packages 

4.1.1 REEDS (Executive Module) 

4.1.2 REEDT 

4.1.3 RCLDS 

4.1.4 RCLDT 

4.1.5 RMTPS 

4.1.6 RMTPT 

4.1.7 RCONS 

4.1.8 RCONT 

4.1.9 RISOS 

4.1.10 RISOT 


SECTION 4.2 

AUXILIARY SUBROUTINE LISTINGS 

4.2.1 

CTLB2 

4.2.2 

MLIB2 

4.2.3 

CAL IX 

4.2.4 

DLINE 

SECTION 4.3 

MAIN MODULE QUESTION FILES 

4.3.1 

7REEDS 

4.3.2 

?RCLDS 

4.3.3 

7RC0NS 

4.3.4 

7RISOS 

SECTION 4.4 

EXAMPLE METEOROLOGICAL DATA F: 
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4.1 


MAIN PROGRAM MODULES AND SUBROUTINE PACKAGES 


This section contains the current listings of current Main Program 
Modules and the respective Main Subroutine Packages. They are listed according 
to the sequences Indicated In Section 4.0. 
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4.2 AUXILIARY SUBROUTINE LISTINGS 

This section contains the current REEDS ai»’ ' 'ary subroutine listings. 
They are listed according to the sequence indicate< . Section 4.0. 
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4.3 MAIN MODULE QUESTION FILES 

This section contdins the Listings of the Main Module Question Files. 
They are listed according to the sequence indicated in Section 4.0 
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4.4 EXAMPLE METEOROLOGICAL SOUNDING FILES 

This section contains the current listings of two (2) example 
meteorological sounding files. They are listed according to the sequence 

Indicated In Section 4.0. 
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